Heisenberg limited measurements with superconducting circuits 
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We describe an assembly of TV superconducting qubits contained in a single-mode cavity. In the 
dispersive regime, the correlation between the cavity field and each qubit results in an effective 
interaction between qubits that can be used to dynamically generate maximally entangled states. 
With only collective manipulations, we show how to create maximally entangled quantum states 
and how to use these states to reach the Heisenberg limit in the determination of the qubit bias 
control parameter (gate charge for charge qubits, external magnètic flux for rf-SQUIDs). 
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The description of the interaction between àtoms and 
quantized modes of the electromagnètic field in a cav- 
ity is called cavity quantum electrodynamics (cQED). 
The first experimental studies used flying Rydberg àtoms 
in a rf resonator With the advent of quantum 

computing several other implcmcntations were developed 
to mimic the quantum properties of àtoms. Among 
those, solid-state implementations are especially interest- 
ing since they offer several advantages over real àtoms: 
these artificial atom properties can be tailored and their 
number and location are fixed. We concentrate our dis- 
cussion on superconducting circuits. Superconducting 
implementations of quantum computing attracted a lot 
of attention in recent years because they are inherently 
sealable and single qubit operations have been demon- 
strated with classical coherent control in a variety of 
qubit types. Several theoretical studies have treated the 
interaction of superconducting qubit (s) with a quantized 
electromagnètic field. The proposal of an on-chip cQED 
experiment using a Cooper-pair box as the artificial atom 
strongly coupled to a one-dimensional cavity .2] is espe- 
cially interesting since it was followed by several experi- 
ments. First, the strength of the coupling was shown to 
be indeed stronger than the different decay constants so 
that the vacuum Rabi splitting was observed Q . Subse- 
quently the ac-Stark shift was measured using a quantum 
non-demolition technique (QND) 4j and the decoherence 
time T2 evaluated from Ramsey fringe experiments . 

In a cavity with a high quality factor, the photons can 
serve as an information bus between several qubits and 
therefore create correlations between distant qubits. Be- 
side its fundamental interest and applications to quan- 
tum information processing, entanglcmcnt offers the ad- 
ditional advantage of allowing improved sensitivity in a 
quantum- limited measurement. Compared to a measure- 
ment made with only one sensor, the sensitivity is im- 
proved by Í/VN when N classical sensors are used to 



measure the same quantity. Now, if N quantum sensors 
are coherently coupled, the sensitivity improvement can 
reach 1/N. This ultimate sensitivity is referred to as the 
Heisenberg limit. Techniques to create and control multi- 
particle entangled states necessary to reach this limit are 
already available in different implementations. For in- 
stanec, interferences between two different polarizations 
(modes) of three and four photon entangled states have 
been observed [Ü, Q and an experiment performed on an 
assembly of three beryllium ions demonstrated a spec- 
troscopic sensitivity improvement Q . Today's solid-state 
sensors do not take advantage of the improvement made 
possible by quantum correlations. Using entanglement in 
a superconducting implementation to enable such Heisen- 
berg limited measurements, could therefore revolutionize 
sensor technology with, for instance, electrometers and 
magnetometers . 

In this work, we use techniques from other fields 
(atòmic QED, ion traps, quantum òptics) to describe 
the interactions and manipulations needed in a super- 
conducting circuit to perform a Heisenberg limited mea- 
surement. We propose to use the photon-qubit inter- 
action to create an effective interaction between distant 
superconducting qubits. This interaction is used to gen- 
erate maximally entangled states which in turn are used 
to beat the Standard quantum limit when measuring the 
natural frequeney of the system. We show how this re- 
sults in a Heisenberg limited estimation of the qubit bias 
control parameter. 

Most superconducting qubits can be described by the 
following single qubit Hamiltonian 
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with the bias B z depending linearly on the dimensionless 
control parameter A: B z = b z (l/2 — X). This Hamilto- 
nian is an approximation vàlid around a symmetry point, 
obtained for À = 1/2, called the degeneracy point. The 
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parameters B Z7 B x , b z and A for the two main types 
of superconducting qubits, can be found in the review 
article on Josephson based devices. The bias is con- 
trolled with an elèctric gate charge n g in a Cooper-pair 
box (CPB) or with an external magnètic flux <j) x in a 
rf-SQUID. 

We now consider that the qubit is contained in a cavity 
described by the Hamiltonian H c = u> c {a?a + 1/2). The 
quantized cavity mode adds an incremental contribution 
SX to the bias so that X + SX = X + A c (a^ + a). The single 
qubit Hamiltonian now reads: 

B z B x &zXc , + \ / \ 

Hq-c = — ~ ~2 Ux + ~2~^ a + a ' Uz ' 

We assume the strong coupling limit g 3> K, 7. We ne- 
glect the cavity decay k and excited state decay 7 for 
the moment. Experiments on a single CPB in a cavity 
support this assumption (g w 10 n). We rewrite this 
Hamiltonian in the eigenbasis of Q which is done by 
performing a rotation of 9 about the y axis. The mixing 
angle 9 is given by tan# = B x / B z . The resulting individ- 
ual contribution is then summed N times to describe an 
assembly of N identical qubits in a single mode cavity: 

H N = ÜS Z + uj c {a)a + 1/2) 
+2 5 (a t + a){S z cos 6 + S x sin9), (3) 

where 2g = -b z X c , Q = ^/B% + B 2 Z and 25 = Yh=i 
Of course, the assumption of exact parameter degener- 
acy is rather stringent. Our intent here is to study the 
ideal case to obtain simple analytic expressions. When 
9 = 7r/2 and the rotating wave approximation is made, 
equation Q takes the form of the A^-particle Jaynes- 
Cummings Hamiltonian. When 9 — 0, equation © de- 
scribes a harmònic oscillator with a conditional displace- 
ment, e.g. a displacement that depends on the total state 
of the system S z . However, the limit 6 1 — > is reached 
when the bias B z is maximum and B x — > 0. In this case, 
the two-level approximation Q , and hence equation , 
is not justified. We assume that the qubit-cavity detun- 
ing A = SI — lü c is much larger than the qubit coupling 
g to the cavity. We write S x = (S+ + S-)/2 and ne- 
glect the rapidly oscillating terms. Hamiltonian can 
be approximately diagonalized with a polaronic transfor- 
mation U — exp ( 9 sl £ 6 (aS+ — aJS—)): 

H N = UH N U^ kÇIS z + Lü c {a)a + 1/2) 

+ X {S 2 -S 2 z +S z + 2a^aS z ). (4) 

We define x — (g sin9) 2 / A. A more complete diago- 
nalization would require to also take into account the 
displaced harmònic oscillator with Ud — exp( 2g ™ s (a — 

af)S z ) and perform the transformation UaU H^U^U d . 
However, this complete transformation introduces terms 
proportional to g 2 /u> c which we neglect because we are 



interested in the regime where g <C A <C uj c . We also 
neglect terms in \ 2 that would appear with the complete 
transformation. Hence, equation J2J is vàlid to first or- 
der in \- The coupling of the qubits with the resonator 
induces a shift in the qubit frequency and a state depen- 
dent shift in the resonator frequency. These Lamb and 
ac-Stark shifts (respectively) were predicted in the case 
of a single CPB in a cavity [2j . 

The novelty of equation J3J compared to previously 
published results is the appearance of an effective inter- 
action H sz = \S 2 between TV qubits mediated by a cav- 
ity photon. This interaction is also known as a one-axis 
twisting interaction because it twists around the z 
axis the quantum fluctuations of the total spin S (for 
a system with several spins). Because of this feature, 
the interaction H sz is a key element in non-optical im- 
plementations of Heisenberg limited estimations. It has 
been used extensively in ion traps exp eriments to gen- 
erate maximally entangled states |lll H2I It has 

been shown to be possible to utilize these states to per- 
form an optimal frequency measurement [Tij and improve 
the estimation of rotation angles 0|. More recently, a 
method involving only collective manipulations has been 
used to perform precision spectroscopy on an assembly of 
six beryllium ions 0] • The first step consists in generat- 
ing the maximally entangled state \ip m ) (see for instance 
reference H3) using the time evolution of H sz over a time 
t sz = ir/2 X : 

IVU =e- i * s ''\-N/2) x = 
±(\-N/2) x + l N + E \+N/2) x ). (5) 

When the number of qubit N is odd, another rotation 
e tf s z j s lieec ied in addition to the e i ^ s ^ 0. We take 
the parity into account in equation (J3J) by setting the 
quantity £ to 2 (1) when N is odd (even resp.). Ini- 
tially all the qubits are in their ground state, so that 
the total wave function describing all the qubits is | J = 
N/2,M = -N/2), = \-N/2) z = \ |) 2 ···| l) N 
(we do not consider the field part of the wave function at 
this point). In order to prepare a quantum state accord- 
ing to equation 1(5} , we need to rotate the initial state 
| — N/2) z around the x axis. Therefore, we define the 
operator Un — e I í Sa: e _lír " í "e~ l ^' Sx for notation conve- 
nience. The average of the total spin vector calculated 
with the wave function \ip m ) is zero, (S) = 0. Thus, the 
natural choice of S z (S) as the observable to be measured 
can not be made. Bollinger et al. showed that the par- 
ity operator YiiLi a zi was an adequate observable that 
could be measured with the state (J5J . However the mea- 
surement of this operator for a large number N of qubits 
is difficult since it requires distinguishing odd and even 
numbers of partides in state | |). A method involving 
only collective manipulations has been proposed in [8| 
to circumvent this problem. First the maximally entan- 
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gled state is constructed. In our case the sequence Un 
is applied to the initial state | — N/2) z . Afterward, the 
system evolves freely during a period T, obeying the dy- 
namics defined by the Hamiltonian Hn- Finally, another 
application of the sequence Un transfers the phase infor- 
mation, N<f>/2, into an amplitudc information of cithcr 
state | + N/2) z or | - N/2) z ; 

\^j) = U N e- lÉNT U N \-N/2) z = 

/N \ 
-ism(j<f>)\-N/2) z + 

i N + E ca S (^<t>)\+N/2) z . (6) 

A measurement will collapse the wave function \ip) on 
either state | + N/2) z or state | — N/2) z with probability 
P T = |(1 + cos(A^)) or P l = 1(1 - cos(A^)). 

We propose to use cavity spectroscopy to infer the 
state of the qubits. Assuming that there is a finite but 
small cavity decay rate k, a signal at frequency uj c sent 
in the cavity will experience a phase shift when it is 
transmittcd. Solving the Heisenberg equation for the 
field creation operator, this phase shift é is given by 
tan $ = ±(2xN) / k. The probability Pf is extracted from 
the time dependence of A measurement scheme as 
been proposed based on this principle to perform a quan- 
tum non-demolition measurement of the state of a single 
Cooper-pair box contained in a cavity ; 2| . The difference 
is that because the coupling to the cavity is y/N stronger, 
the phase shift $ is larger than in the single qubit case. 

The main motivation to use A^-particle maximally en- 
tangled state to perform a spectroscopy measurement is 
to be able to relate the N fold frequency increase to 
the phase uncertainty. The uncertainty on a parameter 
£ can be estimatcd from the error propagation formula 
SÇ = AA/\d(A)/dÇ\ by measuring the operator A. We 
introduce the projection operator A = | + N/2)(+N/2\ 
so that the quantity we propose to measure Pp is the 
average of A over the state of equation JüJ). The vari- 
ance AA 2 is then simply given by Pf(l — Pf) (second 
moment of the Bernoulli distribution) which is equal to 
(s'm(N <fi) / 2) 2 . The denominator of the error propaga- 
tion formula, \dP^/d<p\, is N x sm(Ncj)) /2. Hence, the 
measurement of Pf leads to an estimation of the phase 
uncertainty 5<f> equal to 1/N. The phase acquired during 
the free evolution, by the qubit part of the wave func- 
tion, is (fi + x + 2\n) x T where n is the average pho- 
ton number in the cavity. The contribution of the fre- 
quency shifts is negligible in the expression of the phase 
uncertainty. Thus, the frequency uncertainty is given 
by 5fl = 1/(NT). In superconducting circuits, the pa- 
rameter A controls the level spacing Í2 and therefore the 
uncertainties of both quantities can be related through 
the following relation: 

50 = — y£d=ÍA = 6Jcos6»|5A (7) 
y/W+M 



Hence, a measurement of the frequency íl performed 
with the sequence of operations defined by equation © 
results in a Heisenberg limited measurement of the pa- 
rameter A, e.g. in an improvement of the uncertainty 
SX = 1/(NT x b z \ cos 6*|) associated with the estimation 
of À. This mcthod can be used to improve the estimation 
of the gate charge n g (or the external bias flux cj) x ) in a 
system composcd of N Cooper-pair boxes (or rf-SQUIDs 
resp.) that are coherently coupled. The energies b z and 
B x can be determined from a preliminary spectroscopy 
experiment. Sweeping A when applying a periòdic signal 
will flip the qubits when the frequency is resonant. 

Our scheme is useful away from the degeneracy point 
as at this point cos 8 = 0. However, one should bear in 
mind that the coupling \ decreases as the operating point 
is moved away from the degeneracy point, so a trade-off 
should be made to operate away, while not too far from 
this point. The limited validity range around the degen- 
eracy point is not a limitation in a system composed of 
N Cooper-pair boxes since the different coherence times 
decrease with the distance from this point 0] and there- 
fore the operation far away from it, is not adequate to 
observe coherent cffccts. 

To aid the understanding of our scheme in particular 
and of Heisenberg limited measurements in general, we 
wish to emphasize the difference between the quantum 
observable measured and the parameter that the method 
allows to determino with a better precision. The quan- 
tum observable defines the type of superconducting sen- 
sors or qubit. To simplify the discussion, let's say there 
are mainly two types of superconducting devices, elèctric 
charge and magnètic flux sensors. The quantum observ- 
able measured is then either the elèctric charge n or the 
magnètic flux <fi. Now, the Hamiltonian of the system 
can be tuned or controled with a classical (continuous) 
variable which is the gate charge n g in a charge qubit or 
the external magnètic flux <j> x in a flux qubit. A Heisen- 
berg limited measurement consists of evaluating the un- 
certainty of the parameter (5n g or 5<j) x ) for a given value 
of this parameter with a measurement of the quantum 
observable (n or 4>). The parameters usually estimated 
in a Heisenberg limited measurement are either an energy 
splitting, a rotation angle or a phase delay. Therefore, by 
establishing for the first time a relation between other pa- 
rameters, such as the elèctric charge or the magnètic flux, 
and a quantum measurement, we show how a Heisenberg 
limited measurement can have some applications in sen- 
sor technology. 

In this work, we described a collection of N supercon- 
ducting qubits contained in a single mode cavity. Besides 
the usual shifts in the qubit and resonator freqüències, we 
find that the effective Hamiltonian contains a term x>S 2 
describing the interaction between all the qubits. This in- 
teraction can be used to dynamically prepare maximally 
entangled states. We adapt a method used in ion traps to 
demonstrate the use of these states to reach the Heisen- 
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berg limit in the determination of the qubit frequency. 
Finally, we show that the parameter that controls the en- 
ergy spacing can be estimated with an uncertainty that 
scales inversely with the number N of qubits. Hence our 
work establishes the first formal relation between, either 
the elèctric charge or the magnètic flux in an assembly 
of superconducting devices, and the so-called Heisenberg 
limitcd fluctuations. 
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